Knockdown of Crispld2 in zebrafish identifies a novel network for nonsyndromic cleft lip with or without cleft palate candidate genes.

Orofacial development is a multifaceted process involving tightly regulated genetic signaling networks, that when perturbed, lead to orofacial abnormalities including cleft lip and/or cleft palate. We and others have shown an association between the cysteine-rich secretory protein LCCL domain containing 2 (CRISPLD2) gene and nonsyndromic cleft lip with or without cleft palate (NSCLP). Further, we demonstrated that knockdown of Crispld2 in zebrafish alters neural crest cell migration patterns resulting in abnormal jaw and palate development. In this study, we performed RNA profiling in zebrafish embryos and identified 249 differentially expressed genes following knockdown of Crispld2. In silico pathway analysis identified a network of seven genes previously implicated in orofacial development for which differential expression was validated in three of the seven genes (CASP8, FOS, and MMP2). Single nucleotide variant (SNV) genotyping of these three genes revealed significant associations between NSCLP and FOS/rs1046117 (GRCh38 chr14:g.75746690 T > C, p = 0.0005) in our nonHispanic white (NHW) families and MMP2/rs243836 (GRCh38 chr16:g.55534236 G > A; p = 0.002) in our Hispanic families. Nominal association was found between NSCLP and CASP8/rs3769825 (GRCh38 chr2:g.202111380 C > A; p < 0.007). Overtransmission of MMP2 haplotypes were identified in the Hispanic families (p < 0.002). Significant gene-gene interactions were identified for FOS-MMP2 in the NHW families and for CASP8-FOS in the NHW simplex family subgroup (p < 0.004). Additional in silico analysis revealed a novel gene regulatory network including five of these newly identified and 23 previously reported NSCLP genes. Our results demonstrate that animal models of orofacial clefting can be powerful tools to identify novel candidate genes and gene regulatory networks underlying NSCLP.

Role of WNT10A in failure of tooth development in humans and zebrafish.

BACKGROUND: Oligodontia is a severe form of tooth agenesis characterized by the absence of six or more permanent teeth. Oligodontia has complex etiology and variations in numerous genes have been suggested as causal for the condition. METHODS: We applied whole-exome sequencing (WES) to identify the cause of oligodontia in a 9-year-old girl missing 11 permanent teeth. Protein modeling and functional analysis in zebrafish were also performed to understand the impact of identified variants on the phenotype. RESULTS: We identified a novel compound heterozygous missense mutation in WNT10A (c.637G>A:p.Gly213Ser and c.1070C>T:p.Thr357Ile) as the likely cause of autosomal recessive oligodontia in the child. Affected residues are located in conserved regions and variants are predicted to be highly deleterious for potentially destabilizing the protein fold and inhibiting normal protein function. Functional studies in zebrafish embryos showed that wnt10a is expressed in the craniofacies at critical time points for tooth development, and that perturbations of wnt10a expression impaired normal tooth development and arrested tooth development at 5 days postfertilization (dpf). Furthermore, mRNA expression levels of additional tooth development genes were directly correlated with wnt10a expression; expression of msx1, dlx2b, eda, and axin2 was decreased upon wnt10a knockdown, and increased upon wnt10a overexpression. CONCLUSIONS: Our results reveal a novel compound heterozygous variant in WNT10A as pathogenic for oligodontia, and demonstrate that perturbations of wnt10a expression in zebrafish may directly and/or indirectly affect tooth development recapitulating the agenesis phenotype observed in humans.

Loss of DDRGK1 modulates SOX9 ubiquitination in spondyloepimetaphyseal dysplasia.

Shohat-type spondyloepimetaphyseal dysplasia (SEMD) is a skeletal dysplasia that affects cartilage development. Similar skeletal disorders, such as spondyloepiphyseal dysplasias, are linked to mutations in type II collagen (COL2A1), but the causative gene in SEMD is not known. Here, we have performed whole-exome sequencing to identify a recurrent homozygous c.408+1G>A donor splice site loss-of-function mutation in DDRGK domain containing 1 (DDRGK1) in 4 families affected by SEMD. In zebrafish, ddrgk1 deficiency disrupted craniofacial cartilage development and led to decreased levels of the chondrogenic master transcription factor sox9 and its downstream target, col2a1. Overexpression of sox9 rescued the zebrafish chondrogenic and craniofacial phenotype generated by ddrgk1 knockdown, thus identifying DDRGK1 as a regulator of SOX9. Consistent with these results, Ddrgk1-/- mice displayed delayed limb bud chondrogenic condensation, decreased SOX9 protein expression and Col2a1 transcript levels, and increased apoptosis. Furthermore, we determined that DDRGK1 can directly bind to SOX9 to inhibit its ubiquitination and proteasomal degradation. Taken together, these data indicate that loss of DDRGK1 decreases SOX9 expression and causes a human skeletal dysplasia, identifying a mechanism that regulates chondrogenesis via modulation of SOX9 ubiquitination.

FOXE3 mutations predispose to thoracic aortic aneurysms and dissections.

The ascending thoracic aorta is designed to withstand biomechanical forces from pulsatile blood. Thoracic aortic aneurysms and acute aortic dissections (TAADs) occur as a result of genetically triggered defects in aortic structure and a dysfunctional response to these forces. Here, we describe mutations in the forkhead transcription factor FOXE3 that predispose mutation-bearing individuals to TAAD. We performed exome sequencing of a large family with multiple members with TAADs and identified a rare variant in FOXE3 with an altered amino acid in the DNA-binding domain (p.Asp153His) that segregated with disease in this family. Additional pathogenic FOXE3 variants were identified in unrelated TAAD families. In mice, Foxe3 deficiency reduced smooth muscle cell (SMC) density and impaired SMC differentiation in the ascending aorta. Foxe3 expression was induced in aortic SMCs after transverse aortic constriction, and Foxe3 deficiency increased SMC apoptosis and ascending aortic rupture with increased aortic pressure. These phenotypes were rescued by inhibiting p53 activity, either by administration of a p53 inhibitor (pifithrin-α), or by crossing Foxe3-/- mice with p53-/- mice. Our data demonstrate that FOXE3 mutations lead to a reduced number of aortic SMCs during development and increased SMC apoptosis in the ascending aorta in response to increased biomechanical forces, thus defining an additional molecular pathway that leads to familial thoracic aortic disease.

Lxr regulates lipid metabolic and visual perception pathways during zebrafish development.

The Liver X Receptors (LXRs) play important roles in multiple metabolic pathways, including fatty acid, cholesterol, carbohydrate and energy metabolism. To expand the knowledge of the functions of LXR signaling during embryonic development, we performed a whole-genome microarray analysis of Lxr target genes in zebrafish larvae treated with either one of the synthetic LXR ligands T0901317 or GW3965. Assessment of the biological processes enriched by differentially expressed genes revealed a prime role for Lxr in regulating lipid metabolic processes, similarly to the function of LXR in mammals. In addition, exposure to the Lxr ligands induced changes in expression of genes in the neural retina and lens of the zebrafish eye, including the photoreceptor guanylate cyclase activators and lens gamma crystallins, suggesting a potential novel role for Lxr in modulating the transcription of genes associated with visual function in zebrafish. The regulation of expression of metabolic genes was phenotypically reflected in an increased absorption of yolk in the zebrafish larvae, and changes in the expression of genes involved in visual perception were associated with morphological alterations in the retina and lens of the developing zebrafish eye. The regulation of expression of both lipid metabolic and eye specific genes was sustained in 1 month old fish. The transcriptional networks demonstrated several conserved effects of LXR activation between zebrafish and mammals, and also identified potential novel functions of Lxr, supporting zebrafish as a promising model for investigating the role of Lxr during development.

Regulatory variant in FZD6 gene contributes to nonsyndromic cleft lip and palate in an African-American family.

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect affecting 135,000 newborns worldwide each year. While a multifactorial etiology has been suggested as the cause, despite decades of research, the genetic underpinnings of NSCLP remain largely unexplained. In our previous genome-wide linkage study of a large NSCLP African-American family, we identified a candidate locus at 8q21.3-24.12 (LOD = 2.98). This region contained four genes, Frizzled-6 (FZD6), Matrilin-2 (MATN2), Odd-skipped related 2 (OSR2) and Solute Carrier Family 25, Member 32 (SLC25A32). FZD6 was located under the maximum linkage peak. In this study, we sequenced the coding and noncoding regions of these genes in two affected family members, and identified a rare variant in intron 1 of FZD6 (rs138557689; c.-153 + 432A>C). The variant C allele segregated with NSCLP in this family, through affected and unaffected individuals, and was found in one other NSCLP African-American family. Functional assays showed that this allele creates an allele-specific protein-binding site and decreases promoter activity. We also observed that loss and gain of fzd6 in zebrafish contributes to craniofacial anomalies. FZD6 regulates the WNT signaling pathway, which is involved in craniofacial development, including midfacial formation and upper labial fusion. We hypothesize, therefore, that alteration in FZD6 expression contributes to NSCLP in this family by perturbing the WNT signaling pathway.

Crispld2 is required for neural crest cell migration and cell viability during zebrafish craniofacial development.

The CAP superfamily member, CRISPLD2, has previously been shown to be associated with nonsyndromic cleft lip and palate (NSCLP) in human populations and to be essential for normal craniofacial development in the zebrafish. Additionally, in rodent models, CRISPLD2 has been shown to play a role in normal lung and kidney development. However, the specific role of CRISPLD2 during these developmental processes has yet to be determined. In this study, it was demonstrated that Crispld2 protein localizes to the orofacial region of the zebrafish embryo and knockdown of crispld2 resulted in abnormal migration of neural crest cells (NCCs) during both early and late time points. An increase in cell death after crispld2 knockdown as well as an increase in apoptotic marker genes was also shown. This data suggests that Crispld2 modulates the migration, differentiation, and/or survival of NCCs during early craniofacial development. These results indicate an important role for Crispld2 in NCC migration during craniofacial development and suggests involvement of Crispld2 in cell viability during formation of the orofacies.

MAT2A mutations predispose individuals to thoracic aortic aneurysms.

Up to 20% of individuals who have thoracic aortic aneurysms or acute aortic dissections but who do not have syndromic features have a family history of thoracic aortic disease. Significant genetic heterogeneity is established for this familial condition. Whole-genome linkage analysis and exome sequencing of distant relatives from a large family with autosomal-dominant inheritance of thoracic aortic aneurysms variably associated with the bicuspid aortic valve was used for identification of additional genes predisposing individuals to this condition. A rare variant, c.1031A>C (p.Glu344Ala), was identified in MAT2A, which encodes methionine adenosyltransferase II alpha (MAT IIα). This variant segregated with disease in the family, and Sanger sequencing of DNA from affected probands from unrelated families with thoracic aortic disease identified another MAT2A rare variant, c.1067G>A (p.Arg356His). Evidence that these variants predispose individuals to thoracic aortic aneurysms and dissections includes the following: there is a paucity of rare variants in MAT2A in the population; amino acids Glu344 and Arg356 are conserved from humans to zebrafish; and substitutions of these amino acids in MAT Iα are found in individuals with hypermethioninemia. Structural analysis suggested that p.Glu344Ala and p.Arg356His disrupt MAT IIα enzyme function. Knockdown of mat2aa in zebrafish via morpholino oligomers disrupted cardiovascular development. Co-transfected wild-type human MAT2A mRNA rescued defects of zebrafish cardiovascular development at significantly higher levels than mRNA edited to express either the Glu344 or Arg356 mutants, providing further evidence that the p.Glu344Ala and p.Arg356His substitutions impair MAT IIα function. The data presented here support the conclusion that rare genetic variants in MAT2A predispose individuals to thoracic aortic disease.

Development and plasticity of outer retinal circuitry following genetic removal of horizontal cells.

The present study examined the consequences of eliminating horizontal cells from the outer retina during embryogenesis upon the organization and assembly of the outer plexiform layer (OPL). Retinal horizontal cells exhibit a migration defect in Lim1-conditional knock-out (Lim1-CKO) mice and become mispositioned in the inner retina before birth, redirecting their dendrites into the inner plexiform layer. The resultant (mature) OPL, developing in the absence of horizontal cells, shows a retraction of rod spherules into the outer nuclear layer and a sprouting of rod bipolar cell dendrites to reach ectopic ribbon-protein puncta. Cone pedicles and the dendrites of type 7 cone bipolar cells retain their characteristic stratification and colocalization within the collapsed OPL, although both are atrophic and the spatial distribution of the pedicles is disrupted. Developmental analysis of Lim1-CKO retina reveals that components of the rod and cone pathways initially co-assemble within their normal strata in the OPL, indicating that horizontal cells are not required for the correct targeting of photoreceptor terminals or bipolar cell dendrites. As the rod spherules begin to retract during the second postnatal week, rod bipolar cells initially show no signs of ectopic growth, sprouting only subsequently and continuing to do so well after the eighth postnatal week. These results demonstrate the critical yet distinctive roles for horizontal cells on the rod and cone pathways and highlight a unique and as-yet-unrecognized maintenance function of an inhibitory interneuron that is not required for the initial targeting and co-stratification of other components in the circuit.

Rax is a selector gene for mediobasal hypothalamic cell types.

The brain plays a central role in controlling energy, glucose, and lipid homeostasis, with specialized neurons within nuclei of the mediobasal hypothalamus, namely the arcuate (ARC) and ventromedial (VMH), tasked with proper signal integration. Exactly how the exquisite cytoarchitecture and underlying circuitry becomes established within these nuclei remains largely unknown, in part because hypothalamic developmental programs are just beginning to be elucidated. Here, we demonstrate that the Retina and anterior neural fold homeobox (Rax) gene plays a key role in establishing ARC and VMH nuclei in mice. First, we show that Rax is expressed in ARC and VMH progenitors throughout development, consistent with genetic fate mapping studies demonstrating that Rax+ lineages give rise to VMH neurons. Second, the conditional ablation of Rax in a subset of VMH progenitors using a Shh::Cre driver leads to a fate switch from a VMH neuronal phenotype to a hypothalamic but non-VMH identity, suggesting that Rax is a selector gene for VMH cellular fates. Finally, the broader elimination of Rax throughout ARC/VMH progenitors using Six3::Cre leads to a severe loss of both VMH and ARC cellular phenotypes, demonstrating a role for Rax in both VMH and ARC fate specification. Combined, our study illustrates that Rax is required in ARC/VMH progenitors to specify neuronal phenotypes within this hypothalamic brain region. Rax thus provides a molecular entry point for further study of the ontology and establishment of hypothalamic feeding circuits.

Craniofacial abnormalities result from knock down of nonsyndromic clefting gene, crispld2, in zebrafish.

Nonsyndromic cleft lip and palate (NSCLP), a common birth defect, affects 4,000 newborns in the US each year. Previously, we described an association between CRISPLD2 and NSCLP and showed Crispld2 expression in the murine palate. These results suggested that a perturbation in CRISPLD2 activity affects craniofacial development. Here, we describe crispld2 expression and the phenotypic consequence of its loss of function in zebrafish. crispld2 was expressed at all stages of zebrafish morphogenesis examined and localized to the rostral end by 1-day postfertilization. Morpholino knockdown of crispld2 resulted in significant jaw and palatal abnormalities in a dose-dependent manner. Loss of crispld2 caused aberrant patterning of neural crest cells (NCC) suggesting that crispld2 is necessary for normal NCC formation. Altogether, we show that crispld2 plays a significant role in the development of the zebrafish craniofacies and alteration of normal protein levels disturbs palate and jaw formation. These data provide support for a role of CRISPLD2 in NSCLP.

Genetic ablation of Pals1 in retinal progenitor cells models the retinal pathology of Leber congenital amaurosis.

Mutation of the polarity gene Crumbs homolog 1 (CRB1) is responsible for >10% of Leber congenital amaurosis (LCA) cases worldwide; LCA is characterized by early-onset degenerative retinal dystrophy. The role of CRB1 in LCA8 pathogenesis remains elusive since Crb1 mouse mutants, including a null allele, have failed to mimic the early-onset of LCA, most likely due to functional compensation by closely related genes encoding Crb2 and Crb3. Crb proteins form an evolutionarily conserved, apical polarity complex with the scaffolding protein associated with lin-seven 1 (Pals1), also known as MAGUK p55 subfamily member 5 (MPP5). Pals1 and Crbs are functionally inter-dependent in establishing and maintaining epithelial polarity. Pals1 is a single gene in the mouse and human genomes; therefore, we ablated Pals1 to establish a mouse genetic model mimicking human LCA. In our study, the deletion of Pals1 leads to the disruption of the apical localization of Crb proteins in retinal progenitors and the adult retina, validating their mutual interaction. Remarkably, the Pals1 mutant mouse exhibits the critical features of LCA such as early visual impairment as assessed by electroretinogram, disorganization of lamination and apical junctions and retinal degeneration. Our data uncover the indispensible role of Pals1 in retinal development, likely involving the maintenance of retinal polarity and survival of retinal neurons, thus providing the basis for the pathologic mechanisms of LCA8.

Unique domain appended to vertebrate tRNA synthetase is essential for vascular development.

New domains were progressively added to cytoplasmic aminoacyl transfer RNA (tRNA) synthetases during evolution. One example is the UNE-S domain, appended to seryl-tRNA synthetase (SerRS) in species that developed closed circulatory systems. Here we show using solution and crystal structure analyses and in vitro and in vivo functional studies that UNE-S harbours a robust nuclear localization signal (NLS) directing SerRS to the nucleus where it attenuates vascular endothelial growth factor A expression. We also show that SerRS mutants previously linked to vasculature abnormalities either deleted the NLS or have the NLS sequestered in an alternative conformation. A structure-based second-site mutation, designed to release the sequestered NLS, restored normal vasculature. Thus, the essential function of SerRS in vascular development depends on UNE-S. These results are the first to show an essential role for a tRNA synthetase-associated appended domain at the organism level, and suggest that acquisition of UNE-S has a role in the establishment of the closed circulatory systems of vertebrates.

Role of adenylate cyclase 1 in retinofugal map development.

The development of topographic maps of the sensory periphery is sensitive to the disruption of adenylate cyclase 1 (AC1) signaling. AC1 catalyzes the production of cAMP in a Ca2+/calmodulin-dependent manner, and AC1 mutant mice (AC1−/−) have disordered visual and somatotopic maps. However, the broad expression of AC1 in the brain and the promiscuous nature of cAMP signaling have frustrated attempts to determine the underlying mechanism of AC1-dependent map development. In the mammalian visual system, the initial coarse targeting of retinal ganglion cell (RGC) projections to the superior colliculus (SC) and lateral geniculate nucleus (LGN) is guided by molecular cues, and the subsequent refinement of these crude projections occurs via an activity-dependent process that depends on spontaneous retinal waves. Here, we show that AC1−/− mice have normal retinal waves but disrupted map refinement. We demonstrate that AC1 is required for the emergence of dense and focused termination zones and elimination of inaccurately targeted collaterals at the level of individual retinofugal arbors. Conditional deletion of AC1 in the retina recapitulates map defects, indicating that the locus of map disruptions in the SC and dorsal LGN of AC1−/− mice is presynaptic. Finally, map defects in mice without AC1 and disrupted retinal waves (AC1−/−;ß2−/− double KO mice) are no worse than those in mice lacking only ß2−/−, but loss of AC1 occludes map recovery in ß2−/− mice during the second postnatal week. These results suggest that AC1 in RGC axons mediates the development of retinotopy and eye-specific segregation in the SC and dorsal LGN.

PALS1 is essential for retinal pigment epithelium structure and neural retina stratification.

The membrane-associated palmitoylated protein 5 (MPP5 or PALS1) is thought to organize intracellular PALS1-CRB-MUPP1 protein scaffolds in the retina that are involved in maintenance of photoreceptor-Müller glia cell adhesion. In humans, the Crumbs homolog 1 (CRB1) gene is mutated in progressive types of autosomal recessive retinitis pigmentosa and Leber congenital amaurosis. However, there is no clear genotype-phenotype correlation for CRB1 mutations, which suggests that other components of the CRB complex may influence the severity of retinal disease. Therefore, to understand the physiological role of the Crumbs complex proteins, especially PALS1, we generated and analyzed conditional knockdown mice for Pals1. Small irregularly shaped spots were detected throughout the PALS1 deficient retina by confocal scanning laser ophthalmoscopy and spectral domain optical coherence tomography. The electroretinography a- and b-wave was severely attenuated in the aged mutant retinas, suggesting progressive degeneration of photoreceptors. The histological analysis showed abnormal retinal pigment epithelium structure, ectopic photoreceptor nuclei in the subretinal space, an irregular outer limiting membrane, half rosettes of photoreceptors in the outer plexiform layer, and a thinner photoreceptor synaptic layer suggesting improper photoreceptor cell layering during retinal development. The PALS1 deficient retinas showed reduced levels of Crumbs complex proteins adjacent to adherens junctions, upregulation of glial fibrillary acidic protein indicative of gliosis, and persisting programmed cell death after retinal maturation. The phenotype suggests important functions of PALS1 in the retinal pigment epithelium in addition to the neural retina.

Eye formation in the absence of retina.

Eye development is a complex process that involves the formation of the retina and the lens, collectively called the eyeball, as well as the formation of auxiliary eye structures such as the eyelid, lacrimal gland, cornea and conjunctiva. The developmental requirements for the formation of each individual structure are only partially understood. We have shown previously that the homeobox-containing gene Rx is a key component in eye formation, as retinal structures do not develop and retina-specific gene expression is not observed in Rx-deficient mice. In addition, Rx-/- embryos do not develop any lens structure, despite the fact that Rx is not expressed in the lens. This demonstrates that during normal mammalian development, retina-specific gene expression is necessary for lens formation. In this paper we show that lens formation can be restored in Rx-deficient embryos experimentally, by the elimination of beta-catenin expression in the head surface ectoderm. This suggests that beta-catenin is involved in lens specification either through Wnt signaling or through its function in cell adhesion. In contrast to lens formation, we demonstrate that the development of auxiliary eye structures does not depend on retina-specific gene expression or retinal morphogenesis. These results point to the existence of two separate developmental processes involved in the formation of the eye and its associated structures. One involved in the formation of the eyeball and the second involved in the formation of the auxiliary eye structures.

Regulation and function of foxe3 during early zebrafish development.

In this article, we investigate the expression, regulation, and function of the zebrafish forkhead gene foxe3. In wild type embryos, foxe3 is first expressed in a crescent-shaped area at the anterior end of the prechordal plate, corresponding to the polster. At later stages, the hatching gland, the lens, and the anterior pituitary express this gene. Using morpholinos against the zinc finger Kruppel-like factor 4 (KLF4) we show that foxe3 is regulated differently in the polster and in the lens. In the absence of KLF4, expression of foxe3 in the polster is not activated, whereas in the lens placode the expression of KLF4 is not required for the transcription of foxe3. The expression of foxe3 is also regulated by the hedgehog and nodal signaling pathways. foxe3 expression is altered in the hedgehog pathway mutants iguana and you-too and the nodal pathway mutant cyclops. foxe3 function is necessary for the execution of lens-specific gene expression and lens morphogenesis, as the knockdown of foxe3 results in a loss of platelet-derived growth factor receptor alpha (pdgfralpha) expression and in the vacuolization of the lens.

Rx-Cre, a tool for inactivation of gene expression in the developing retina.

Rx is a homeobox-containing gene that is critical for vertebrate eye development. Its expression domain delineates a field of cells from which the retina and the ventral hypothalamus develop. The 5' upstream regulatory sequences of the medaka fish Rx gene are functionally conserved during evolution to a degree that they direct gene expression into the Rx-expressing field of cells in mice. Using these sequences, we made a Cre line that can be used for inactivation of gene expression in the developing retina.

Identification of 315 genes essential for early zebrafish development.

We completed a large insertional mutagenesis screen in zebrafish to identify genes essential for embryonic and early larval development. We isolated 525 mutants, representing lesions in approximately 390 different genes, and we cloned the majority of these. Here we describe 315 mutants and the corresponding genes. Our data suggest that there are roughly 1,400 embryonic-essential genes in the fish. Thus, we have mutations in approximately 25% of these genes and have cloned approximately 22% of them. Re-screens of our collection to identify mutants with specific developmental defects suggest that approximately 50 genes are essential for the development of some individual organs or cell types. Seventy-two percent of the embryonic-essential fish genes have homologues in yeast, 93% have homologues in invertebrates (fly or worm), and 99% have homologues in human. Yeast and worm orthologues of genes that are essential for early zebrafish development have a strong tendency to be essential for viability in yeast and for embryonic development in the worm. Thus, the trait of being a genetically essential gene is conserved in evolution. This mutant collection should be a valuable resource for diverse studies of cell and developmental biology.

Expression of FoxE4 and Rx visualizes the timing and dynamics of critical processes taking place during initial stages of vertebrate eye development.

Several transcription factors have a critical function during initial stages of vertebrate eye formation. In this paper, we discuss the role of the Rx subfamily of homeobox-containing genes in retinal development, and the role of the Foxe3 and FoxE4 subfamily of forkhead box-containing genes in lens development. Rx genes are expressed in the initial stages of retinal development and they play a critical role in eye formation. Elimination of Rx function in mice results in lack of eye formation. Abnormal eye development observed in the mouse mutation eyeless (ey1), the medakatemperature-sensitive mutation eyeless (el), and the zebrafish mutation chokh are caused by abnormal regulation or function of Rx genes. In humans, a mutation in Rx leads to anophthalmia. In contrast, Foxe3 and FoxE4 genes are expressed in the lens and they play an essential role in its formation. Mutations in the Foxe3 gene are the cause of the mouse mutation dysgenetic lens (dyl) and in humans, mutation in FOXE3 leads to anterior segment dysgenesis and cataracts. Since Rx and FoxE4 are expressed in the earliest stages of retina and lens development, their expression visualizes the timing and dynamics of the crucial processes that comprise eye formation. In this paper we present a model of eye development based on the expression pattern of these two genes.